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We propose a new process of /i~e~ — > e~e~ in a muonic atom for a quest of charged lepton flavor 
violation. The Coulomb attraction from the nucleus in a heavy muonic atom leads to significant 
enhancement in its rate, compared to /i + e~ — > e + e~. The upper limit of the branching ratio is 
estimated to be of the orders of O(10 -17 - 10 -18 ) for the photonic and the four Fermi interactions 
from the present experimental constraints. The search for this process could serve complementarily 
with the other relevant processes to shed lights upon the nature of charged lepton flavor violation. 
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Charged lepton flavor violation (cLFV) is known to 
be one of the important rare processes to search for new 
physics beyond the Standard Model (SM). Various theo- 
retical models predict sizable rates of cLFV processes, 
which are just below the present experimental upper 
limits. The on-going and future experiments for cLFV 
searches would reach sensitivities in the range of pre- 
dictions by many theoretical models. At this moment, 
the cLFV searches with muons present the best limits 
owing to a large number of muons available for measure- 
ments [H- Typical cLFV processes with muons include 
fi + — > e + 7, /i + — > e + c + c~ and pT - e~ conversion in 
a muonic atom (/-i~N — > e~N). However, after the dis- 
covery of cLFV process in future, many other different 
cLFV processes should be studied to shed lights upon un- 
derstanding of the nature of the cLFV interactions and 
develop insights into new physics responsible for cLFV. 

In this letter, we would like to propose a new cLFV 
reaction process of a bound pT in a muonic atom, which 
is 

p~e~ — » e~e~, (1) 

where p~ and e~ in the initial state of Eq.(JT|) are the 
muon and the atomic IS electron(s) bound in a Coulomb 
field of the nucleus in a muonic atom respectively. 

This p~cT — » e~e _ process in a muonic atom has 
various significant advantages. First of all, this pro- 
cess could have not only photonic dipole interaction but 
also four-Fermi contact interaction, as in the processes 
of p + — > e + e~e~ and — > e~N, whereas p + — > 0+7 
has only the former. This would potentially allows us to 
investigate the full structure of new physics beyond the 
SM. Secondly, this process has a two-body final state, in 
which the two signal electrons are emitted almost back 
to back and each of them has an energy of about a half 
the muon mass, m M /2. This would provide a cleaner ex- 
perimental signature as well as a larger final-state phase 
space than p + — > e + e + e~ decay. Also, in comparison 
with the p + — » c + 7 search, the measurement of this pro- 
cess would be relatively easier since no photon detection 



is involved. Thirdly, one can consider a similar reaction 
process with a muonium, such as p + e~ — > e + e _ . How- 
ever, the rate of this /i + c~ — » c + e~ process can not be 
large because of small overlap between the p + and e - 
wave functions. However, in a muonic atom of atomic 
number Z, we can increase the overlap between the p~ 
and e _ wave functions if an atom of large Z is chosen. 
The enhancement occurs owing to the Coulomb interac- 
tion from the nucleus which attracts the IS state electron 
wave function towards the p~ and the nucleus. The ex- 
pected rate would increase by a factor of (Z — l) 3 . For 
example, the rate for a lead [Z = 82) is 5 x 10 5 times 
that of the /i + e~ — > c + c~ reaction. However, in a muonic 
atom, a nuclear muon capture process occurs in addition 
to the normal Michel muon decay. But since a lifetime 
of a muonic atom changes from 2.2 ps to ~ 80 ns from 
a hydrogen to a lead, the reduction of the number of 
muons of a factor of at most 20 can only be expected. 
Therefore, the net increase of the branching ratio would 
become significant for a large Z . A potential disadvan- 
tage is that the rate of the reaction process like this might 
not be large enough compared to rare cLFV muon decays. 
Therefore, in this letter we will evaluate the rate of this 
p~e~ — > e~e~ process, and discuss its related issues. 

We describe the process of p~e~ — > e~e~ in a muonic 
atom of Eq. ([1]) by an effective Lagrangian at the energy 
scale of the muon mass m M . Following Ref . [l[ , we define 

^-e-^e-e- = [«V MR ^"^-^ 

+ m M A L JlZa^ e K F^ 

+ 9i(WeL)(ene L ) + 32(757e R )(er;eR) (2) 
+ 33(MR7' 1 eR)(eR"7 A1 e R ) + g 4 ('pL'y t "ei J )(eZj fl ei J ) 
+ 95 (mrT^ e R ) (^l7/j e L ) + i g6(7iL7 A ' e L)(eR7MeR) 
+ (H.c.)] , 

where Gf — 1.166 x 10~ 5 GcV~ 2 is the Fermi coupling 
constant, and Al, r and gj's (i £ {1,2,- •• ,6}) are di- 
mensionless coupling constants. The first two terms in 
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FIG. 1: The process /i~e~ — ► e~e~ induced from the photonic 
interactions. The black dot indicates the effective interaction 
that is absent from the Standard Model. 



the brackets of Eq. ^ are the photonic interaction con- 
tributing to the process of Eq. (|TJ) through the diagrams 
shown in Fig. [TJ The remaining terms are for the direct 
four-Fermi contact interaction. 

We estimate the branching ratio for the process of 
Eq. ([T]), which takes place in a muonic atom with an 
atomic number Z. The initial state has the muon and 
the electron in their IS ground state of the atomic orbits. 
For simplicity, we ignore the three-momenta of the initial 
bound muon and electron. The final state has the two 
electrons, which, in the lowest order, can be treated as 
monochromatic plane waves that propagate with oppo- 
site momentum vectors. Each of the two electrons in the 
final state takes energy of about m M /2, when the bound 
effects at the IS state and the Coulomb interaction from 
the nucleus can be neglected. 

We begin with the first case where the four-Fermi 
interaction is dominant and Ar and are negligibly 
small, making no contribution of the photonic interac- 
tions. The remaining four-Fermi interaction allows the 
processes such as /j.~e~ — > e~e~ and /i + — > e + e + e~. The 
cross section of the process of Eq. ([TJ is calculated to be 
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where G = G\i 

Gij = \gi\ 2 + \gj\ 2 and G ■ ■ = Re(g*gj) . The transition 
rate is then given by 
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Here we took into account the facts that the IS state 
can accommodate two electrons, and that the nuclear 
charge is shielded by the negative muon. We used the 
non-relativistic wave functions given by 
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(5) 



where r is the radial coordinate, so that ip^J{0; % ~ 1) = 
[(Z - l)am c ] 3 / 2 /y/n . The rate of Eq. @ is enhanced for 



a larger atomic number, Z, by a factor of (Z — l) 3 , giving 
a notable advantage for heavy nuclei. This enhancement 
comes from the factor of '(0; Z — 1)| 2 , and the large 
positive charge of a heavy nucleus strongly attracts the IS 
wave functions of the lcptons toward the nucleus position, 
rendering the overlap of the two wave functions large, and 
enhances the transition of the process of Eq. ([Ij. We 
normalize the rate of Eq. (QJ by the lifetime of a muonic 
atom, fu, to define the branching ratio of this process as 
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The value of f M ranges from 7^ = 2.19x10 6 s for 1 H to 
r p = (7 - 8) x l(T 8 s for 238 U as listed in Ref. 0. This 
is shorter than the lifetime of free muons, r M = 2.197 x 
10" 6 s [3], which is equal to 1927r 3 /(G|m^) at the lowest 
order. The obtained branching ratio of Eq. (jSJ) is to be 
compared with that of //+ — > e + e + e~, the process that 
arises from the same elementary process. This branching 
ratio with no photonic interaction are given by Q 
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; (G 12 + 16G 3 4 + 8G 56 ). (7) 



The contribution from the interference among the four- 
Fermi interactions are not present in Eq. ([7]) , whereas it is 
found to be present in Eq. ^ as the terms of G^-'s. The 
absence of the interferences is due to the large momenta 
of the final electrons. The search for the process of Eq. (UJ) 
will thereby serve complementarily with that for /i + — > 
e + e + e~. The ratios of the two branching ratios of 

Br(u~c~ — >• c~c~) „ , .„ Q/TO r \ 3 f„ , . 
5 , . + _ <192rr (Z - 1 3 a 3 (^ -± , 8 
Br(/i+ — >• c+c+e ) ^m^J 

is obtained by assuming G/(G\2 + I6G34 + 8G56) ~ O(l)- 
Hence, the upper limit on Br(/i~e~ — > e~e~) is con- 
strained by the existing limit of Br(/x + — > e + e + e~) < 
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Figure [5] shows these upper limits as a function of an 
atomic number Z by the dotted curves. These upper 
limits are obtained by taking the present experimental 
limit of S max = 1.0 x 10" 12 from the SINDRUM ex- 
periment The light-shaded region in Fig. [5] is ex- 
cluded by the current SINDRUM limit. The reciprocal 
of the shown limit gives an estimation of the number 
of muons that is required to detect events of the pro- 
cess of Eq. (UJ). Let us take an example of the gold 
atom (Z = 79): our estimation requires a collection of 
(4.21 x lO" 19 )- 1 = 2.38 x 10 18 muon events to improve 
the current limit of S max = 1.0 x 10~ 12 . Assuming the 
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FIG. 2: Limits on the branching ratio to the process of Eq. |T]) 
imposed by the limits on Br(/i+ — » e + 7) and on Br(/i+ — > 
e+e+e - ). Models without photonic interactions are excluded 
in the light-shaded region. Similarly, models with photonic 
interactions are excluded in the dark-shaded region by the 
current experimental limits. 



detection efficiency of O(10%), the required number of 
muons amounts to a few times 10 19 . Compared to this 
number is O(10 18 -10 19 ) of muons, which is the goal of 
the highly intense muon beamsplanned in the near fu- 
ture in search for the cLFV [6j-|8|]. We thereby find that 
the current limit could be reachable within the capabil- 
ity of these muon sources. Note that relativistic effects 
and the binding energy of the leptons may be consequen- 
tial for heavy nuclei, but are not included in our present 
estimation. 

Let us now turn to the second case where the pho- 
tonic interaction is present and it dominates over the 
four-Fermi interactions. In this case, the cross section, 
rate, and branching ratio of the process of Eq. ((TJ) are 
calculated to be 
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and 



Br(/i c — » e e ) 
= 1536^ 2 (Z - 1)V(|A R | 2 + L4 L | 2 )^^ . 



(12) 



respectively. On the other hand, the branching ratio of 
fi + — > c+c+c~ is given by Q 
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We then have 

Br(/i~e~ — » e - c~) 
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which we combine with the upper limit Br(/i+ — > 
e+e+e") < B max to yield 
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This upper limit is overlaid in Fig. [5] by a dash-dotted 
curve, according to the aforementioned SINDRUM limit 
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The presence of the photonic interactions gives rise to 
another cLFV process /i+ — » c+7 as well, and search for 
this process also put a limit to Br(/i~e~ — > e~e~). The 
branching ratio to this process is given by 
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which we compare with Eq. (1121) as 
Br(/i - e~ — > e~e~ 



Br(/i+ -> e+7) 



4 (Z-l) 3 a 4 ^^ 
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Then the limit on Br(/i c — > e e ) is estimated from 
Br(/x+ e+7) < B max as 
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A solid curve in Fig. [5] presents the upper limits given 
in Ea.(fl8]) with B max = 1.2 x 10~ n , which is the cur- 
rent upper limit from the MEGA experiment. A dashed 
curve in Figure [5] also shows the limits with the B max — 
1.7 x 10 -13 , which is the goal value of the MEG exper- 
iment [L0| |. Even the current limit on the branching ra- 
tio to /i+ — > e+7 overwhelms the limits on Br(/Lt _ e _ — > 
c~c~) when the photonic interaction is dominant. Ac- 
cordingly, the present excluded region is above the cur- 
rent MEGA limit. Let us estimate, as we did earlier, the 
required number of muons to detect events of the process 
of Eq. ([TJ), taking an example of the gold atom (Z = 79): 
an accumulation of (1.044 x IO -17 ) -1 = 9.58 x 10 16 muon 
events is necessary to surpass the sensitivity of MEGA, 
and (1.480 x IO" 19 )- 1 = 6.76 x 10 18 to exceed that of the 
MEG goal. The required number of muons are estimated 
to be 10 18 and a few of 10 19 , respectively, assuming the 
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0(10%) of the detection efficiency again. These are not 
capable now but would be possible in future with the 
planned highly intense muon sources. 

The expected magnitudes of the branching ratio of the 
process of Eq. (fT|), which is driven by both the four- 
Fermi interaction and the photonic interaction, are found 
to be not significantly large even with the enhancement 
of (Z — l) 3 . Thus, this would not be the first process 
to cultivate the discovery frontier of the cLFV searches. 
However, thanks to the enhancement of (Z — 1) , this 
process can be accessible in future by next-generation 
high-intensity muon beams to produce muons of an order 
of O(10 18 -10 19 ) per year. On the other hand, the search 



for it 



would not be possible with these in- 



tensities. Even higher beam intensities of O(10 21 ) muons 
per year are envisioned in the future accelerator projects 
such as muon colliders, neutrino factories, and Higgs fac- 
tories 



A number of critical issues should be carefully stud- 
ied to consider experimental feasibility of this search. 
First, one should take account of potential modifica- 
tion of energy spectra of each of the emitted signal elec- 
trons due to the bound effects in a muonic atom. It 
has been known that the electron spectrum from bound 
muon decays in a muonic atom is strongly deformed 



for a large atomic number Z [12j. For the process of 
Eq. fl}, however, it is speculated that the sum of the 
two signal electrons are still fixed to the initial energy of 
?n M + m c — (binding energies) due to the small kinetic en- 
ergy of the recoiled nucleus. Secondly, the wave functions 
of the initial and final leptons should be improved to take 
account of the nuclear charge distribution and relativistic 
effects. A binding energy of the initial muon need to be 
also taken into account. Thirdly, the interactions of the 
initial leptons need detailed treatments. Various differ- 
ent models may predict the dependence of the rate on the 
initial spin state. And the interaction between the two IS 
electrons in the atomic orbital of the muonic atom can be 
looked at. Possible contributions from 2S, 3S, • • • states 
should also be considered. Finally, an important experi- 
mental issue is the estimation of backgrounds, including 
the physical SM backgrounds from yTe~ — > e~e~vu de- 
cay, and accidental backgrounds that are known to be 
detector-dependent. We will discuss these issues in our 
future works. 

In summary, the new cLFV process /i~e~ — > c~c~ 
in a muonic atom is proposed. This process has the 
rate enhancement of (Z — l) 3 over the fx + e~ — » e + e~ 
owing to the Coulomb interaction from the nucleus in 
a muonic atom. This process has a final state of two 
electrons, which would be experimentally very clear sig- 
nature. The upper limits of the branching ratio of the 
orders of 0(1CP 17 -10~ 18 ) are estimated separately for 
the photonic and the four Fermi interactions from the 
other cLFV experimental results. Once this process is 



observed, CP violation might be studied by comparing 
this process with /i + — > e + e + e~. 
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